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ABSTRACT 

We provide a novel method to infer finger flexing motions using a four-channel 
surface electromyogram (EMG). Surface EMG signals can be recorded from the 
human body noninvasive^ and easily. Surface EMG signals in this study were 
obtained from four channel electrodes placed around the forearm. The motions consist 
of the flexion of five single fingers (thumb, index finger, middle finger, ring finger, and 
little finger) and three multi-finger motions. The maximum likelihood estimation was 
used to infer the finger motions. Experimental results have shown that this method can 
successfully infer the finger flexing motions. The average accuracy was as high as 
97.75%. In addition, we examined the influence of inference accuracies with the 
various arm postures. 
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INTRODUCTION 

Compared to the invasive type electromyogram 
(EMG) for measuring each motor unit action po- 
tential (MUAP) on muscle fiber in a small area, 
surface EMG signals can be recorded from the 
human body non-invasively and easily over wide 
skin areas. To extract the motion information, 
surface EMG signal processing has been studied, 
especially focusing on controlling electronic devices 
such as a human supporting prosthetic arms and 
smart interface (Sornmo and Laguna, 2005). Since 
motor commands are transmitted through internal 
nerve tissues or muscles, inferring such commands 
only from the surface EMG signals is not an easy 
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task. 

Thus many attempts to extract motor information 
from surface EMG have been made. A study of 
Englehart et al. (2001) used wavelet analysis and 
PCA analysis to classify two arm motions and two 
wrist motions. Also, Englehart and Hudgins (2003) 
used the absolute mean value, the zero crossing 
rate, and the wavelength to classify four arm and 
wrist motions. Momen et al. (2007) used RMS and 
FCMs to classify arm and wrist motions. Hudgins et 
al. (1993) obtained an ensemble average of EMG 
signals, and used a neural networks (NN) to 
classify four arm motions. 

While the surface EMG-based motion inference 
has been focused on upper limb or hand motions, 
studies on inferring dexterous finger motions have 
not been done so much. Nishikawa et al. (1999) 
used the Gabor transform and absolute mean value 
to extract the features and classify four wrist 
motions and six finger motions in real-time learning 
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Fig. 1. Electrodes orientation for recording EMG on the forearm. 



based on NN. Nagata et al. (2007) used absolute 
sum analysis, canonical component analysis, and 
minimum Euclidean distance to classify four wrist 
motions and five finger motions. Uchida et al. 
(1992) used FFT analysis and NN to classify four 
finger motions. Chen et al. (2007) used mean 
absolute values (MAV), the ratio of the MAVs, 
autoregressive (AR) model, and linear Bayesian 
classifier to classify 5~ 16 finger motions. However, 
NN classifiers (Hudgins et al., 1993; Chen et al., 
2007; Nagata et al., 2007) have strong dependency 
on initial parameter conditions. Also, linear classi- 
fiers (Uchida et al., 1992; Englehart et al., 2001; 
Englehart and Hudgins, 2003; Momen et al., 2007; 
Nagata et al., 2007) are simple but show low 
accuracy compared with nonlinear classifiers. 

In this paper, we propose a novel method to infer 
the finger motions using surface EMG signals. We 
recorded the EMG signals during finger flexing 
motions. To infer the finger motions, we used the 
maximum likelihood method. The likelihood was 
obtained by building probabilistic models of EMG 
activities. The information entropy of the EMG mag- 
nitude was calculated to quantify the EMG acti- 
vities. Experimentally we demonstrated that the 
proposed method could infer the finger flexing mo- 
tions including multi-fingers as high as 97.75%. 

In addition, we examined the influence of 
inference accuracies with various arm postures. 

Multi-channel EMG recording 

A signal acquisition device (Poly-G-A, Laxtha Inc., 
Korea), bi polar Ag-AgCI electrodes (Dual electrode 
#272, Noraxon U.S.A. Inc.), and snap electrodes 
(SECS-4, Laxtha Inc., Korea) were used. The EMG 
signal was observed at the sampling frequency of 
512 Hz and was filtered to reduce power line noise. 

Two able-bodied subjects without any upper limb 
deficiencies participated in performing motions for 
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Fig. 2. Relaxation and finger flexing motions (Fi~F 8 ). 

this study and was trained before signal recording. 
The subjects were instructed to maintain their own 
flexing speed and strength. 

The EMG signal was acquired with four channels 
on the left forearm (around the wrist) for 3.0 se- 
conds. Fig. 1 shows electrode orientation for EMG 
recording. During the first 1.5 seconds, the subject 
did not perform any motion - termed 'relaxation' - 
and waited for an acoustic alarm. After 1.5 se- 
conds, the subject heard the acoustic alarm and 
then flexed a finger immediately. After 3.0 seconds, 
the EMG signal recording procedure ended. In the 
experiment, the subject performed eight finger flex- 
ing motions including three multi-finger motions. Fig. 
2 shows the detail of the finger motions. The 
subject took a long and sufficient rest between 
recording procedures to avoid muscle fatigue or 
cramp. 

PROPOSED METHOD 

Preprocessing of EMG signal 

Here we consider four EMG channels and eight 
finger flexing motions. Let c and k be the channel 
and motion indexes, respectively. Then during F k 

motion, the recorded raw EMG from the C th chan- 
nel can be denoted by r k [n] where n is a discrete 
time index. 

Since the EMG signal has both positive and 
negative levels, we take an absolute of r k [n\. 
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Fig. 3. Examples of raw EMG signal and absolute of EMG signal 
corresponding the 2 nd channel during the F 3 motion. The re- 
corded samples that treated as signals are limited in 1.5-2.2 
sec (unshaded area). 



1 

H{X)= Yi p(x m )log 2 —7- } 

m = l P\ x m' 



(2) 



where X denotes discrete random variable, 
p(x m ) means the probability of the random variable 
X when X equals x m and s" =l p{x m ) = l. For 
c h channel and F k motion, we define the pro- 
bability p k (m) as 



P k (m)-- 



# of samples £= /„, 
Pr of total samples 



/ =i x [•]!- m ~ l < x [-}<X — 1 



m = 1,...,M, 



(3) 



\r h n 



(1) 



For example, Fig. 3 shows the surface EMG 
signal of the 2 nd channel during the F 3 motions, 

r 3 [n] and x 3 [n]. 



Statistical EMG modeling 

The amplitude of the EMG signal is a funda- 
mental quantity which increases monotonically with 
the force developed in the muscle (Sornmo and 
Laguna, 2005). Moreover, the EMG signal repre- 
sents a stochastic signal (Sornmo and Laguna, 
2005). 

The entropy (information entropy, H) was used to 
represent the characteristics of each EMG signal 
activity during finger flexing motions. The entropy is 
a measure that can reveal how much dynamical 
nature and the information (Shannon and Weaver, 
1963) which the signal contains. Some EEG an- 
alysis methods based on entropy have been de- 
veloped. (Capurro et al., 1999; Martin et al., 2000; 
Tong et al., 2002; Shin et al. 2006; Sabeti et al., 
2009) Although the amplitude range or the spectral 
bandwidth of EMG and EEG signal are remarkably 
different, but they have a non- stationary stochastic 
property in common. 

Conventionally the entropy is calculated as 
follows: 



where M determines the number of total bins 
and x max indicates the largest magnitude of the 
signal acquisition device. This value is set as 
x max = i,050yuV in the study but can be adjusted 
to correspond to the subject's EMG characteristics. 
Then, the information entropy of the EMG signal is 
given by 



m = l Pky i > 



(4) 



After accumulating the entropy in eq. (4), we 
build the probability density function of the entropy 
based on the Gaussian model: 



fk(h) = 



1 



Ittct l 



-exp 



1 



lot 



(5) 



»k = e I h k] < 



Here f k (h) is the probabilistic model based on 
the entropy of EMG signals. The example of this 
modeling is given in Fig. 4. Then we can get the 
probability density functions f k on each channel for 
8 different finger motions. Fig. 5 shows the resulting 
models of 8 finger motions and 4 channels. As- 
suming the statistical independence in each chan- 
nel's EMG signal, the likelihood function L{k) be- 
comes 
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Fig. 4. Building the probabilistic model of the entropy of the EMG signals. 
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Fig. 5. Probabilistic model of entropy 
the EMG signals. 



L{k)=f[f k {h). (6) 

c=l 

Then, we find k which maximizes L(k) such that 

k = arg max L(k) (7) 

k 

where k means the estimated motion. Fig. 6 
shows the likelihood for various candidates of finger 

motions. 



EXPERIMENTAL RESULTS 

A cross-validation is used to estimate the in- 
ference accuracies. The number of training data is 
25 and that of the test data is 25. All test data are 
excluded from the training data. The results are 
averaged over 500 combinations. 

The accuracy for six motions (F v F 2 , F 3 , F 4 , 
and F 8 ) using the proposed method was 100.00% 
without any confusions as shown in Table 1. The 
proposed method can infer the motions more 
accurately than the studies of Uchida et al. (1992) 
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Fig. 6. Likelihood for various can- 
didates of finger motions (each value 
is normalized so that the maximum 
was equal to one. k and k denote 
the actual and the inferred motions, 
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Table 1. Inference accuracy (%) for six motions (subject A, fc: 
actual motion, fc: inferred motion) 



Table 2. Inference accuracy (subject A, test data: 25, training 
data: 25) 
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Fig. 7. 

25). 
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Confusion matrix (subject A, test data: 25, training data: 



and Chen et al. (2007) 

Fig. 6 shows Likelihood for various candidates of 
finger motions (each value is normalized so that the 

maximum was equal to one. k and k denote the 
actual and the inferred motions, respectively). 



Table 3. False table (subject A, fc: actual motion, fe: inferred 
motion) 



0.0% 
2.0±1.4% 



8.0±2.7% 



5.0±2.2% 



3.0±1.7% 
0.0% 
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(e) k = F„ k=F 5 (f) k=F„ k=F 6 

(g) fc = F 7 , k= F 7 (h) fc = fc= F 8 

The average inference accuracy for four motions 
(flexion of thumb, flexion of index finger, flexion of 
middle finger, flexion of all fingers) in a study 
(Uchida et al., 1992) was 86% and the accuracy 
for six motions (extension of thumb, extension of 
index finger, extension of middle finger, extension 
of ring finger, extension of little finger, and 'hook' 
gesture) was 96.83% (Chen et al., 2007). 

The inference results for eight motions (F 1 ~F S ) 
are summarized in Table 2 showing the high 
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Fig. 8. Defined basic various nine postures. 
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Fig. 9. Inference accuracies on various postures, Pi~P 9 (subject B). 



inference accuracy. The inference accuracies on 
F v Fi, F 5 , and F s are 100% but the slight 
decreases on F 2 , F 3 , F 6 , and F 7 arise from false 
inference. Fig. 7 shows the confusion matrix which 



visualizes the false inference among the performed 
motions and Table 3 shows the false inference rate 
which indicates the rate when performed motion is 
inferred incorrectly. The false inference rate on F 2 
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Table 4. Average Inference accuracy 

Training dataset 

Test dataset 





Pi 


P2 


P 3 


P 4 


P 5 


Pe 


P7 


Ps 


P 9 


Pi 


74.72% 


26.83% 


27.12% 


37.07% 


30.45% 


38.34% 


27.49% 


20.83% 


40.83% 


p 2 


38.24% 


83.78% 


41 .78% 


49.20% 


62.29% 


43.65% 


43.20% 


63.49% 


62.94% 


Ps 


53.79% 


44.59% 




64.21 % 


43.02% 


42.75% 


46.85% 


48.74% 


54.43% 


P4 


54.69% 


50.54% 


61 .09% 


81.15% 


50.16% 


56.67% 


40.70% 


45.14% 


46.19% 


p 5 


39.58% 


65.70% 


41 .63% 


49.42% 


83.60% 


44.83% 


42.45% 


47.22% 


34.72% 


p 6 


39.38% 


37.42% 


40.77% 
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40.42% 


62.62% 


35.79% 


38.02% 


36.75% 
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41 .20% 


41.14% 


32.22% 


38.51 % 


42.24% 


40.30% 




43.98% 


39.84% 


p 8 


43.42% 


65.82% 


40.42% 


35.68% 


52.72% 


37.87% 


44.91 % 


76.83% 


50.05% 


p 9 


45.33% 


50.74% 


36.34% 


40.93% 


44.00% 


34.34% 


37.70% 


44.58% 


77.55% 



confused with F 6 is 2.00±1.41%. The false in- 
ference rate on F 3 confused with F 6 is 5.00±2.20% 
while that on F 6 confused with F 3 is 8.00±2.74%. 
The false inference rate on F7 confused with F 3 is 
3.00+1.72%. 

Influence of human arm posture 

We set nine arm postures to find out how much 
accuracies hold down when human changes arm 
posture. Fig. 8 shows basic postures. 

Matrices on the diagonal site in Fig. 9 show 
inference accuracies for five single-finger flexion 
motions of each posture; the number of training 
data was 25 and that of test data was equal. Each 
experiment was done 50 times. All test data were 
excluded from the training data. The average 
inference accuracies are 74.72%, 83.78%, 76.22%, 
81.15%, 83.60%, 62.62%, 63.07%, 76.83%, and 
77.55%, respectively. 

We tested again with the data of another pos- 
tures as a training dataset. Other matrices in Fig. 9 
show the results. The average inference accuracies 
are shown in Table 4. The shaded cells indicate the 
maximum average accuracy in the same row and 
these match the diagonal cells. 

Accordingly, the data of an posture is not ap- 
propriate to the reference of other postures. This 
fact indicates that a robust finger motion decoding 
method cannot be implemented without a scheme 
detects changes of arm's posture. 

DISCUSSION 

We have proposed a new method to infer finger 
flexing motions based on the entropy and the 



maximum likelihood estimation. The average recog- 
nition accuracy for six finger motions was 100.00% 
and the accuracy for eight finger motions was 
97.75%. Also we have shown the quantitative 
interpretation of the need of avoidance of use of 
reference dataset cross arm posture because that 
influences the correlations between the EMG sig- 
nals and the performed finger flexing motions. This 
study may trigger a new type of human-computer 
interface with user's intuitive hand motions which 
could control electronic devices. 
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